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Abstract 

We study the first sub- leading correction 0((lns) ) to the cusp anomalous dimension in the 
high spin expansion of finite twist operators in M = 4 SYM theory. Since this approximation is 
still governed by a linear integral equation (derived already from the Bethe Ansatz equations 
in a previous paper), we finalise it better in order to study the weak and strong coupling 
regimes. In fact, we emphasise how easily the weak coupling expansion can be obtained, 
confirms the known four loop result and predicts the higher orders. Eventually, we pay 
particular attention to the strong coupling regime showing agreement and predictions in 
comparison with string expansion; speculations on the 'universal' part (upon subtracting the 
collinear anomalous dimension) are brought forward. 
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1 Plan of the work 



Let us consider a (massless) gauge theory with a certain field in the adjoint representation (or 
q in the fundamental one). To fix ideas, we may focus our attention on the planar sl(2) sector of 
TV = 4 SYM with local operators 

Tr(^0 L ) + ...., (1.1) 

where T> is the (symmetrised, traceless) covariant derivative acting in all possible ways on the L 
bosonic fields <fi. At finite twist L the high spin behaviour of the anomalous dimension 

l(g,s,L) = A(g;s,L)-L-s, (1.2) 

will in general be determined by the universal (in the sense here of being twist and flavour inde- 
pendent [U though theory dependent) scaling function f(g) (cf. for instance 0, HI Q], EJ El [2] 
and references therein): 

y(g, s, L) = f(g) In s + f sl (g, L) + . . . , (1.3) 

(with a factor 1/2 for the field q in the fundamental representation [5]) with this parametrisation 
of the 't Hooft coupling A = 87r 2 g 2 . The universal scaling f(g) equals in the L = 2 sector twice 
the cusp anomalous dimension of the light-like Wilson loops [I] and therefore inherits this name. 
We denote above by dots terms which are going to zero, whilst we wish to investigate the sub- 
leading (constant) scaling f s i(g,L) (for QCD see also [7]). In fact, the leading scaling f(g) has 
been constrained (for all the values of g) in [2] (cf. also [5] for the relation with the Bethe root 
density) by a linear integral equation derived from the asymptotic Bethe Ansatz [H] and from the 
same Bethe Ansatz the sub- leading one f s i{g, L) has been shown in [10], after careful consideration 
pTj . to enjoy a linear integral equation with the same kernel (but different inhomogeneous term; 
for an alternative approach using the non-linear integral equation [T2], see [H]). It is, of course, 
non-universal in the sense again that it depends explicitly on the twist L (and, if there are other 
representations than the adjoint one, on the representation carried by the field, as well as on the 
theory). Yet, it will be conjectured in the following with some evidence to be wrapping-free (at 
least for the smallest twist L = 2), thus sharing the destiny of f(g) in this respect q The main 
evidence is coming from the agreement with the recent string computations [H], which should 
seriously allow for gauge loop wrapping effects, as they live at very large coupling g. In addition, a 
very recent analysis [T3] a la Luscher has proved for the scalar twist two operator that the wrapping 
correction at four loops goes to zero (like In 2 s/s 2 ). 

In any case, thanks to Kotikov, Lipatov, Onishchenko and Velizhanin [16] it is by now well 
established that (in Af = 4 SYM) 

fM = f.i(g,2) (1.4) 

is the same for the operators with the smallest twist L = 2 (and determines also that of some twist 
three operators, see p2]). in fact, all twist two operators belong to the same supermultiplet, and 
their anomalous dimension is expressed in terms of a 'universal' function defined by the scalars 

1 As already clear, but also stressed in the following it is not universal as f(g). 
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luniv.(s) = 7^(s) and with shifted Lorentz spin for gauginos 7^(s) = 7um«( s + 1) an d gauge fields 
7 A (s) = 7«n™(s + 2). 

Still at the leading twist L = 2 Dixon, Magnea and Sterman [18] have pointed out that, 
although not universal (as f(g)), its first logarithmic integral, X \g), equals the sub-leading 
correction to the logarithm of the (diverging) scattering amplitude (the following is just a sketchy 
formula whose precise definitions and regular isat ion rely on [18] or references therein, in particular 
for the connexion between fg^(g) and f s i(g) §) 

ln ^/^_4!M + ...., (1 . 5) 

e 2 e 

in dimensional regularisation to D — 4 — 2e, plus the first logarithmic integral of a universal function 
h(g) (times a representation depending factor we omit here because not extant in SYM4). In terms 
of the by-logarithm- derivative function of f^, (<?), i.e. the so-called collinear anomalous dimension 
fsi'(g)) this means the important relation between sub- leading term 

fsi(g) = fsi'(g) + h(g) . (1.6) 

In other words, this [L = 2) sub-leading scaling in the anomalous dimension (11.31) contains a 
universal part, h(g), which, once subtracted from it, yields (upon integration) the sub-leading cor- 
rection in the amplitude (jl.5P ; both the sub-leading terms are not universal, albeit their difference 
ought to be. As we will see later (cf. (I3.12p ). the linear integral equation leads us to a split form of 
f s i(g,L) into a twist-dependent part and another independent of it, f s i(g) (cf. (13.121) below). The 
latter contains somehow the universal h(g) and thus might share with it, at very strong coupling, 
its logarithmic behaviour (cf. below (15. 2p and (15.41) which are in agreement with the string theory 
calculation by [19]). This may lead us to the possible proposal 

h(g) = f(g)[\n-^--( 1 -- 3 -\n2) + ...}, (1.7) 

with an intriguing object between square parenthesis (and the presence of f(g) may be perhaps 
motivated by the relation of h(g) to a Wilson loop [H]). In fact, at the first order, recalling that 
f(g) = 2\p2g + . . . , we would obtain 

h(g) = 2y/2g [In -L - (I - I l n2 ) + . . . ] , (1.8) 

which coincides with the tested function in string theory for gluons and quarks [20J. 

Note added: An interesting paper [2T] appears today in the web archives. It seems to contain 
an analysis of the strong coupling expansions of f s i(g) in overlap with some of the results of this 
work. 



2 The same connexion holds between /^^(g) and f s i'(g)- A similar one between /( 2 \g) and f(g) (cf. 
references therein). 
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2 One loop results 



The present analysis fully relies upon the main achievements of [TO] , where it was shown that the 
relevant integral equations which enter the computation of the anomalous dimension at large spin 
are linear. Let us briefly remind some results which will be useful later. 

For what concerns the one loop results, we consider equation (3.52) of that paper. The meaning 
of such equation is that, in the high spin limit s — > oo, the one loop density of roots is approximated, 
up to orders o(s°), by the function whose Fourier transform &o(k) reads as 



f (l - e~^J + e"^ M - cos 



(j (fc) = -2tt^ '- ^ ^ - ATr5(k) In 2 + o(s°) . (2.1) 

sinh 

We also showed that the one loop energy, 



/°° dk 
— e(k)& (k) — (L — 2)e(0) 



with the function 

e(u) = — T e(k) = 2ne~^ , (2.2) 

in the high spin limit behaves as 

E (s, L) = 4 In s - 4(L - 2) In 2 + 4 7s + o(s°) , (2.3) 
and the one loop anomalous dimension reads 

7o (a, L) = g 2 E (s, L) = 4g 2 In s - [4(L - 2) In 2 - 4 lE ]g 2 + g 2 o(s°) . (2.4) 

3 All loops results 

In order to study the higher than one loop density of roots <J^(fc) it is convenient to define the 
quantity 

2 sinh ^ 

S(k) = ^^-a B (k) . (3.1) 
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We then consider equation (4.11) of |10) 3 I . Passing to Fourier transforms, we obtain that the 
function S(k) satisfies the linear equation 

S(k) = A[i _ M V2 gk)] + _L £j ^ [f>(-ir+i J r (V2 g A;) J r (V2y/i) 1 " ^^ e-l + 

oo oo 

+ sgn(/i) ^ ^+1+2,(5) (-l) r+1/ e"^ (jr^iVZgtyJr+iviV^gh) - J r _i (V2gh) J r+2u (V2gk)j ■ 



r=2 u=Q 



S , (/i)-47rln2 5(/i)-7r(L-2)- ^- - 2tt ■ + o(s°) , (3.3) 



sinh ^ sinh ^ sinh ^ 

where the functions c r ^ s {g) enter the definition of the "dressing factor" (for their definition, see 
e.g. [2]). As follows from the generalisation of the Kotikov-Lipatov identity [8], the relation of the 
function S{k) with the (all-loops) energy (anomalous dimension) is 

<y(g,s,L) = 2S(0). (3.4) 

We now want to put equation (13.31) in a form which is suitable for analysis of both the weak and 
the strong coupling limit. It is convenient to restrict to the domain k > and expand S(k) in 
series of Bessel functions 

S(k) = £sM*yM (3.5) 

P =i 

with the energy H3.4[) given by 

1 (g,s,L) = V2gS 1 (g). (3.6) 

On the other hand, for what concerns the coefficients S p (g), after some simple calculations we find 
the following system of equations 

S 2p ^(g) = 2V2g(\n s + lE )5 p , + 4(2p - 1) J™ f ^-i (^) _ 

- 2(2p-l)(L-2) r dh J 2p MV2gh) _ ^ Z^^S^g) , 

Jo h t^i 

(3.7) 

S 2p (g) = 4 + 8p/ f^M + 2(L-2)-4,(L-2) f + 

he— 1 j ft e2 + l 



^2p,2m-l (g)S 2m {g) 



m=l m=l 
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The functions ^0(^)1 Fh(u), appearing in (4.11) of [10j . are related to <Jo(u) and <jh{u) by the relations 

«7 (u) = -j-FoM , a H (u) = -rF H (u) . (3.2) 
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In writing (13. 7p we used the notations 

J ( , T ( \ x x y m H dh J n (V2gh)J m (V2gh) 

J2p-i(x) = J2p~i{x) - 6 Pil - , Z n>m (g) = J — e h _i • ^ 3 - 8 ) 

The structure of the forcing terms appearing in the linear system (13.71) suggests that its solution 
and - consequently - the (all-loops) anomalous dimension can be split as 

S r (g) = S* ES {g) Ins + (L - 2) SV>(g) + S e r xtra (g) , (3.9) 

where S EES (g) is the (coefficient of the) solution of the BES equation [2], Sr (g) is the (coefficient 
of the) first generalised scaling density [221 1221 I2U [23 I2E] and the extra coefficients S^ xtra (g) are 
solutions of the following system 

S$Z(g) = 2vW P ,! + 4(2 P - 1) jf f - 

oo 

- 2(2p-l) J2 Z2 P -i,m(9)S e ^ ra (g) , 

771=1 

(3.10) 

sir (9) = 4+8 p/°°f^f + 

oo oo 

+ a p Z2 P ,2m-i{g)s^-M -4pY1 Z2 P ,2 m {g)SlT{g) ■ 

777=1 777=1 

From the splitting (13. 9p of their solution, we obtain that the all loops energy at high spin is 

j(g : 8l L) = f(g) In s + (L — 2)f^(g) + f sl (g) + o(s°) , (3.11) 

where f(g) is the cusp anomalous dimension, f^(g) is the first generalised scaling function and 
fsi(g) = V2gSl xtra (g) comes from the solution of (I3.10p . In other terms, the sub- leading correction 
is worth 

f sl (g,L) = f sl (g) + (L-2)fW(g), (3.12) 

with explicit mention to its twist dependence (since f s i(g) does not depend on L). 

In the next sections we will analyse the weak and the strong coupling expansion of f s i(g, L). 



4 Weak coupling 

Since the system ( 13. 10p is linear (as the original linear integral equation [10]), it can be systemati- 
cally and very easily expanded at small g, giving the weak coupling (convergent) series of f s i(g, L) 
(via (I3.12p ). This kind of computations can be automatised by using a program for the symbolic 
manipulation in order to go on to an arbitrary order @ . 

4 In what follows we report the results up to six loops, but we easily reached loop ten. We invite the interested 
reader to contact the authors for having these expressions, since their writing would fill some pages in. 
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In view of (13.121) . it is natural to start by the twist two formula up to the order g 12 : 



4 

9 



+ ^( 7 r 2 C(3) + 30C(5)) (-|) 6 + 

(7A(3) + 50vr 2 C(5) + 2625C(7)) (^L 
128 R _ , 832 . 560 



+ 



( — A(3) + 192C(3) J + -^W) , 



-tt 2 C(7) + 14112C(9)) (^) " + 



+32 (164C(3) 2 C(5) + 4851C(H)) ) +.... (4.1) 

And we end up with the general twist expansion (we also expanded f^(g) according to [131 [10] ) 



i 



fsi(9, L) = f sl (g) + (L - 2)f^(g) = ( lE — (L — 2) \n2)f{g) + 8(2L - 7)C(3) I ^= ) 

~l {* 2 am - 4) + 3(2LL - 62)C(5)) " + 

8 / \ 8 

+ - (tt 4 C(3)(3L - 13) + 75(46L - 127)C(7) + 5(111 - 32)vr 2 C(5)) UL) + 

i^7r 6 C(3)(HL - 49) + 8 (2695C(9)L + 16C(3) 3 L - 7154((9) - 56C(3) 3 ) + 
+ ^(25L - 64)vr 2 C(7) + A(i 03 L - 310)vr 4 C(5) J " + 

(32 8 
-vr 4 C(7)(295L - 772) + -tt 2 (l519C(9)L + 24C(3) 3 L - 3626C(9) - 88C(3) 3 ) 

+8 (33285C(11)£ + 536C(3) 2 C(5)L - 85974C(11) - 1728C(3) 2 C(5)) + 

8 , , „ , , 8(2956L - 13231)tt 8 C(3)\ / g 

+ 2023L - 6266 A 5 + — ; sv ; 

945 V ; SV ; 14175 / V\/2 
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(4.2) 



Importantly, this expansion for twist L = 2, 3 agrees with those up to four loops present in the 
literature [TBI l2Tl [T3"] and checked by different means (cf., for instance, [Hj for a summary). Thanks 
to [151 EE] we can also argue that for twist two (at least up to four loops) f s i(g) enjoys the property 
of being wrapping free, and we would like to conjecture the same in general about fsAg, L). 

In fact, last but not least, after the completion of this work we became aware lj that a five 



3 We ought to thank Matteo Beccaria for this private communication. 



7 



loop calculation (wrapping inclusive) on twist L = 3 would exactly confirm the expansion of this 
section. 



5 Strong coupling 

Let us pass to consider the strong coupling limit of (13.101) . From the results in Appendix A, we 
get that the leading contribution at large g in the forcing terms of (I3.10P is 

-2y/2g\ng6 Ptl , (5.1) 

which allows to write down the following analytic leading behaviour 

fM = -f(9)hig + .... (5.2) 

where f(g) is the cusp anomalous dimension at strong coupling [2"§l [3"0] and the dots denotes 
contributions of order g or smaller. This important point will be clarified in the remaining part of 
the section. 

In order to analyse the next to leading terms of f s i{g) at large g, we solve numerically the linear 
system (13.101) . Indeed, this system is amenable for the usual numerical treatment as developed in 
[31J. The kernel is the same as in the study of the universal scaling function, settling any issue 
about convergence. What changes now is the nature of the forcing terms, which is expected to 
produce some subtle effects at strong coupling. In particular, as explained in appendix A, the 
forcing term which contains J\(x) shows a logarithmic leading asymptotic behaviour for g — > oo. 

A (numerically) clean way to study the system at large g is to explicitly subtract such a 
logarithmic dependence in equations (I3.10p . We are then led to consider the objects S^ xtra (g), 
satisfying the linear system 



S-Jto) = 2V~2 9 - In *f) S rA + 4(2p - 1) [ 

OO 

- 2(2p-l)J2Z2 P -i,m(9)S^ ra (g), 



m=l 



(5.3) 



sir (9) = 4 + 8p/ ^My^) + 



m=l m=l 

With such a subtraction we have a simple expression for f s i(g), 



f sl (g) = f(g) In ^ + V2 g S^g) , (5.4) 
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where f(g) is the usual cusp anomalous dimension and Sl xtra (g), supposed to be free of logarithms, 
can be easily studied numerically (a plot of the complete function f s i(g) in the range g G [0, 5] can 
be found in fig. 1 together with a magnification of the weak-coupling behaviour). From a general 
point of view, it is likely to expect that \f2g S{ xtra (g) can be expanded in powers of g as follows, 

V2gSr tra (g) = k l9 + k + — + 0(1/ g 2 ), g - oo . (5.5) 

9 



The usual numerical analysis [31] . including up to L = 30 Bessel functions in the Neumann 
expansion, gives the following best fit coefficients: 

k num = _ 2 .828426 ± 0.000001; k% um = 0.3238 ± 0.0001; k™\ m = -0.01194 ± 0.00015. (5.6) 

Inserting this expansion into (15.41) we obtain 



f sl (g) = 2V2g 



2V2 31n2 , 2V2 c K 2^2 k-i , „,ln<? 
In a -= — In 1 = In 1 =— + U{ 

9 2V2ng 9 2V2ng 8n 2 g 2 g 2y/2g 2 9 
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(5-7) 



where K = (3(2) is the Catalan's constant and c\ = —k 1 /(2\/2), cq = koir. Now, we may also infer 
the plausible analytic values^ 

4K-9(ln2) 2 , . 

Cl = l, co = 6ln2-vr, k-x = 7=^ ^ = -0.0118253.... (5.8) 

4y 27r 2 

In particular, for twist two this value of cq = 6 In 2 — 71 entails the non-published evaluation 
c = 6 In 2 + 7r by N. Gromov, reported in the added note of |14j . 

Eventually, for generic twist L plugging into (13.121) the asymptotic value of the first generalised 
scaling function '(g) = — 1 + 0(e ^5) [22], we end up with 



f sl (g,L) = 2V2g 



, 2^ 3 In 2 , 2^2 Co + (2 — L)ix K 2v^2 lnp, 
In ci -= — In 1 -= In 1 + 0(- 



g 2V2ng g 2y/2irg 8n 2 g 2 g 2^2 g 2 g 3 

(5.9) 

This shows clearly how, in this asymptotic expansion, the only trace of the twist comes up in 
the piece Co + (2 — L)ir and thus cancels out completely, at order O(s ), in the asymptotic (large 
g) expansion of A — s = 7 + L. This cancellation is due to the simple asymptotic expansion 
f (9) = ~ l + 0(e~^) [22] and holds for any L. Furthermore, it makes possible that the constant 
term in A — s = 7 + Lat order O(s ) be - = 6 ln . 2+7r for any twist and thus future comparison 
with string theory (which does not distinguish between null and small values of L |14j). 



6 Summarising 

In the present work we have performed a detailed study of the sub-leading (constant) contribution 
to the anomalous dimension of the twist operator in the large spin expansion. 

6 They seem also to be in agreement with the expansion by Frcyhult and Zieme |21j 
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Figure 1: Plot of f s i(g) with magnification for g e (0, 1). 

To this aim, we fully exploited the results obtained in [TU] in order to write down a linear 
integral equation for f s i(g,L). This approach allowed us to study its weak-coupling expansion up 
to the order g 20 (and even further), and to access the strong coupling behaviour by means of a 
modification of well established numerical techniques. An entire range plot is depicted by Fig. 1. 

The strong coupling regime is of particular interest, because it shows a leading logarithmic 
behaviour (with the 't Hooft coupling), that we are able to single out analytically by inspecting 
the structure of the inhomogeneous term in the linear integral equation (see appendix A). Moreover, 
relying on the numerical analysis we are in the position to disentangle the terms of the asymptotic 
expansion; in particular, we confirm the string result in [T3] , additionally fix the constant c for the 
twist 2 case and understand how this constant is the same for any other twist. 

It is also important to stress that our computations, based on the asymptotic Bethe Ansatz, 
are in full agreement with previous findings [161 123 021 EH HS1 [28] . 

Finally, our formalism also allows the treatment and can be seen as a particular case of the 
high spin and large twist case, when the ratio between the logarithm of the spin and the twist is 
fixed: this is the subject of the contemporaneous paper [32J. 
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A Strong coupling of the extra terms 

We briefly analyse the strong coupling (asymptotic) behaviour of the forcing terms of the system 
(13.101) for the extra coefficients S^ xtra (g). From the asymptotic formula 



00 







dhh . J J^M = ± 2m+ ;-'^ , r+s>1 , (a.i) 



where B n are the Bernoulli numbers, we deduce the large g behaviour of all the forcing terms in 
(13.101) . but the one involving J\. The strong coupling limit of such term can be evaluated after 
writing it as 

dx J 2 (x) + J (x) - 1 
2 e 75^ _ 1 

Then we compute separately the strong coupling limits of the addends containing J 2 and Jo — 1, 
respectively. For the former, we use (lA.ljl . For the latter, one uses the integral representation 

J (x) - 1 = — / d6 (e lxsine - 1) (A.3) 

and then integrates first on x, getting an expression in terms of ip functions, and finally on 9, after 
developing the ip functions for large argument. With this procedure we get the asymptotic formula 

[~dh MV2gh) g g g (\ \ 1 fl\ 

We see that such term at large g behaves as — ^ ln^. Therefore, it dominates the strong coupling 
behaviour of the forcing terms of (I3.10P . 
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